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ABSTRACT 

We present a simultaneous analysis of 18 galaxy lenses with time delay measurements. For 
each lens we derive mass maps using pixelated simultaneous modeling with shared Hubble con- 
stant. We estimate the Hubble constant to be 66^4 km s _1 Mpc _1 (for a flat Universe with 
Vtm = 0.3, Qa — 0.7). We have also selected a subsample of five relatively isolated early type 
galaxies and by simultaneous modeling with an additional constraint on isothermality of their 
mass profiles we get Ho = 76^ km s _1 Mpc _1 . 



Subject headings: gravitational lensing, cosmological parameters 



1. Introduction 

The Hubble constant is one of the most impor- 
tant parameters in cosmology. It determines the 
age of universe, the physical distances to objects, 
constrains the dark energy equation-of-state and 
furthermore, it is used as a prior in many cosmo- 
logical analysis. Hence, it is essential f or cosmol- 
ogy to know the precise value of Ho (|Riess et al.l 
2009). The newest measurements give a quite well 
defined Ho with estimated errors at the 5-10% 
level but unfortunately Ho differs among different 
cosmological methods and there is only marginal 
consistency between their la errors (see Table 1). 

Moreover, those methods are based on differ- 
ent physical principles and more importantly they 
measure different consequences of Ho. Both SNe 
and Cepheids measure luminosity distance but at 
different scales (distant and local Universe respec- 
tively). Cepheids can provide a luminosity dis- 
tance via the period-luminosity relation but only 
in the local Universe, on the other hand, the signif- 
icantly brighter SNe with their characteristic peak 
luminosity can measure cosmological distances but 
need to be calibrated with Cepheids. The SZ ef- 
fect, which is based on high energy electrons in 
a galaxy cluster distorting the CMB through in- 
verse Compton scattering, is proportional to the 
gas density of the galaxy cluster and that com- 
bined with the cluster's X-ray flux gives an esti- 



mate of the angular diameter distance. Finally, 
the angular power spectrum of the CMB gives 
information about many cosmological parameters 
which have correlated effects on the power spec- 
trum, especially strong degeneracies are between 
Ho and Q m and £V Furthermore, there are sig- 
nificant differences in the obtained results even 



within one method (e g . SNe iBranch et al.l ([ 19961 ): 
Sandage et all (|2QQ6h : iRiess et al.l ([2009l Tsee Ta- 



ble 1). 

Additionally, any astrophysical approach suf- 
fers from systematic uncertainties, e.g., super- 
novae are not standard but standardized can- 
dles with possible redshift evolution, CMB have 
various parameter degeneracies and interference 
with foregrounds, and the SZ method is assum- 
ing spherical symmetry on often significantly non- 
spherical clusters of galaxies. It is therefore impor- 
tant to explore complementary methods for mea- 
suring Ho- Gravitationally lensed quasars QSOs 
offer such an attractive alternative. 

As shown by iRefsdall (|l964h the Hubble con- 
stant can be measured based on the time delay 
At between multiply lensed images of QSOs, be- 
cause Ho oc 1/At, provided that the mass distri- 
bution of the lens is known. Time delays measure 
D °Dls° S ' wnere ^ol, Ads, Als are the angular 
diameter distances between observer and lens, ob- 
server and source, and lens and source, respec- 
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tively. Gravitational lensing has its degeneracies 
but it is based on well understood physics and 
unlike distance ladder methods there are no cali- 
brati on issues ([Branch et al.|[l996[ ISandage et al.l 

2qq3). 



Gravitational lensing has, up to now, not been 
seen as reliable as other leading cosmological 
methods. Determination of the Hubble constant 
using lensing is problematic, because the mass dis- 
tribution of a lens strongly influences the result of 
Ho and unfortunately we never have a complete 
knowledge of that, hence a choice of a lens model 
is needed. 

Recently, however, time-delay lenses have suc- 
cessfu lly been used for Ho estimation. In partic- 
ular, Oguril (j2007f ) used a Monte Carlo method 
to combine lenses and derive d Hp = 70 ± 6 
km s _1 Mpc _1 . ISaha et al.l ([20061 ) obtained 
Hq = 72in km s _1 Mpc _1 using a combina- 



tion of 10 lenses, IColesI (|2008l ) got H = 71± 6 8 
km s ~ 1 Mpc~ 1 using a combination of 11 lenses, 



69-7Z5.5 km 



and ISuvu et al.l (j20Q9h found H - « Q 7+4,9 
s _1 Mpc _1 by detailed analysis of one gravitational 
lens, B1608+656. 



This p aper extends the work of IS aha et al 



([20061 ) and lColesI ([20081 ). where results on Ho were 
presented using combined modeling of 10 and 11 
lenses, respectively. We have used those systems 
with refined properties of lenses as a part of our 
sample and have added new systems that have 
been discovered and monitored during the past 4 
years. We therefore now posses an almost dou- 
bled sample of systems with measured time delay, 
which demonstrates that, gravitational lensing is 
valuable method for Ho estimation. 



Table 1: Ho comparison 



Method 



Author 



CMB (3 years) 
CMB (5 years) 
la SNa 
la SNa 

SZ 
Cepheids 



3.2 
+2.6 



73.2 

71.9_ 27 

74.2 ±3.6 

62.3 ± 1.3 

65t? 
72 ±8 



Soergel et al. (2007) 
Komatsu et al. (.2009) 

Riess et al. ^2009) 
Sandage et al. (2006) 

Jones et al. (2005) 
Freedman et al. (2001) 



2. Pixelated modeling 

Two different approaches for modeling lenses 
are commonly used. The first one, the non- 
parametric (pixelated) method, generates a large 
number of models which perfectly fit the data, 
each of them giving a different result which can 
then be averaged. For the second method, model 
fitting, one assumes parametrized models of the 
mass distribution of the lens. 

Pixelated modeling has the advantage of allow- 
ing the lens shape and profile to vary freely. It 
does not presume any parameters and can provide 
models that would not be possible to reproduce 
with parametric modeling. Using this approach 
to a combined analysis of a large sample of lenses 
is a powerful solution to the modeling problem in 
gravitational lensing. 

In this p aper we use the no n -analy tical method 
created by ISaha Williams! (|2004f ) - PixeLens. 
PixeLens generates an ensemble of lens models 
that fit the input data. Each model consists 
of a set of discrete mass points, the position of 
the source and optionally, if the time delays are 
known, Ho. The time delay At is the combined 
effect of the difference in length of the optical path 
between two images and the gravitational time di- 
lation of two light rays passing through different 
parts of the lens potential well, 



Ar = 



1 + z L DqsDql 



9(6) 



Here and j3 are the positions of the images and 
the source respectively, z\, is the lens redshift, and 
9 is the effective gravitational potential of the 
lens. 

The arrival time at position is defined in Pix- 
eLens, as a modeled surface, 



r(0) = l -\o\ 



\n\Q-O f \K0')d 2 O f , (2) 



where k is surface mass density. 

The errors of the positions of observed images 
and redshifts of source and lens are of the order 
of a few percent, thus can be ignored. The main 
source of errors in the data comes from time delays 
between images. 

We have used PixeLens to generate a set of 100 
models for a sample of lensing systems (see §3). 
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PixeLens produces an ensemble of models with 
varying Ho, each consisting of sets of mass pixels, 
which exactly reproduce the input data. More- 
over, it also models several lenses simultaneously 
enforcing shared Ho for all lenses. The images po- 
sitions, the source and lens redshifts, and the time 
delay are assumed to be accurate enou gh for their 
errors to be ignored (jSaha et al.l 120061 ). PixeLens 
also imposes secondary constraints on mass maps: 
non-negative density; smoothness, where the den- 
sity of a pixel must be no more than twice the av- 
erage density of its neighbors; the mass profile is 
required to have 180° rotation symmetry (except 
if it appears very asymmetric); the shear is allowed 
within 45° of the chosen direction; circularly av- 
eraged mass profile should nowhere be shallower 
nor steeper than r~ amin and r~ amax respectively, 
where a m i n and a max is defined in PixeLens as the 
minimum and maximum steepness, those values 
can be chosen by the user however, the default Pix- 
eLens constraint is minimal steepness a m i n = 0.5; 
and finally, additional lenses as point masses can 
be constrained. PixeLens does not use flux ra- 
tios as constraints becau se of the possible influ - 
ence of redden ing by dust (lElfasdottir et al . 2006), 
microlensing ( Paraficz et al.l 20061) or small-scale 
struc ture in the lens potential ([Dalai fc Kochanek 
2002h . 



3. Data set 

To date, there are 19 gravitational lens sys- 
tems with published time delays. Table 2 sum- 
marizes the information about these 19 systems. 
We have made an attempt to use all the conclu- 
sions previously drawn about their shape, external 
shear, profile, etc. We have used the newest/best 
measurement of positions and redshifts of images 
and lens. Apart from the main lensing galaxies 
we have also included all the galaxies that might 
contribute to the lensing. We added them when- 
ever they are visible in the field. These systems 
are RX J0911+055, HE 1104-181, SBS 1520+530, 
B1600+434 and B1608+656. All the mass maps 
of the doubly imaged quasars are required to have 
180° rotation symmetry and in case of the quadru- 
ply imaged systems we allow the lens to be asym- 
metric if it has been reported asymmetric, which 
is the case for: HE 0435-1223, SDSS J1004+4112, 
RX J1131-1231, B1608+656. A constant external 
shear is allowed, for the lenses where the morphol- 



ogy shows evidence of external shear or the exis- 
tence of external shear has been reported, which is 
the case for HE 0435-1223, RX J0911+055, FBQ 
J0957+561, SDSS J1004+4112, PG 1115+080, 
RX J1131-1231, SDSS J1206+4332, B1600+434, 
SDSS J1650+4251, WFI J2033-4723. 

One lens systems has been excluded from 
our a nalysis, B1422+231. iRavchaudhurv et al. 
(|2003l ) indicated that the time-delay m e asure - 



ments made by Patnaik fc N arasim hal (|2Q0lh 
are p ossibly inaccurate. IPatnaik fc Narasimhal 
(l200lh reported ATi 2 = 7 .6 ± 2.5 days, whereas 



Ravchaudhurv et al.l (j2QQ3h lens m odeling predicts 
AT12 = 0.4/z -1 days. According to IRavchaudhurv et al, 
(2003) this value would not be expected to show 
up in the IPatnaik fc Narasimhal (j200lf ) data, 
which sampled ever y 4 d ays. We follow the 
Ravchaudhurv et "all (|2QQ3h prediction also be- 



cause in our analysis the system gives an unrea- 
sonably low Hubble constant Ho — 12 + 3 km 
s _1 Mpc _1 , and we exclude the system in what 
follows. 

Figure 1 shows a mosaic of the average mass 
distributions for the remaining 18 lenses. 

4. Full set results 

Our resulting Ho distribution is shown in Fig- 
ure El We have performed the calculation of 18 
systems for a flat Universe with Q\ = 0.7, Q m = 
0.3 and we obtain Ho = 66^4 km s _1 Mpc _1 at 
68% confidence and Ho = 66^7 km s _1 Mpc _1 at 
90% confidence. We also note that for (zl) = 0.6, 
(zs) = 1.8, that are average lens and source red- 
shifts of our sample, the inferred Ho should in- 
crease by 2% for an open Universe with Q\ = 0.0, 
Q m = 0.3 and decrease by 7% for an Einstein-de 
Sitter Universe = 0.0, Vt m = 1.0. 

Figure 3 presents the comparison between es- 
timation from our sample of lenses and samples 
from previous lensing studies using 15 lense£| 



Excluding B1422+231. 
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Table 2 

Properties of time-delay lenses 



System 






P.A. 


Point mass a 






Ar c 


Reference 












(") 


(") 


(days) 








PG 1115+080 


0.311 


1.722 


-45° 





0.381 


1.344 




1,2, 


3,4 












-0.947 


-0.69 


13.315;* 
















-1.096 


-0.232 



















0.722 


-0.617 


11.711;* 






HE 2149-275 


0.495 


2.030 


_ 





0.714 


-1.150 




5, 6, 7 












-0.176 


0.296 


103 ± 12 






SDSS J1206+4332 


0.748 


1.789 


50° 





0.663 


-1.749 




8, 


9 












0.566 


1.146 


116±| 






HE 0435-1223 


0.454 


1.689 


-65° 





-1.165 


0.573 




10, 


11 












1.311 


-0.03 


9 -.+0.78 

A- -L _q y\ 
















0.311 


1.126 


6 +l."07 
— 1 .08 
















-0.226 


-1.041 


Q Q7+1.31 

— 1.67 






SDSS J1650+4251 


0.577 


1.547 


80° 


o 


0.017 


0.872 




12, 


13 












-0.206 


-0.291 


49.5±1.9 






QSO J0957+561 


0.356 


1.410 


-30° 





1.408 


5.034 




14, 15, 16 












0.182 


-1.018 


422.6 ± 0.6 
















2.860 


3.470 



















-1.540 


-0.050 









RX J1131-1231 


0.295 


0.658 


-35° 





-1.984 


0.578 




17 


,3 












-1.426 


-1.73 


i r+2.49 
2 02 
















-2.016 


-0.610 


9.61+H? 

— 1.57 
















1.096 


0.274 


87 + 8 






RX J0911+055 


0.769 


2.800 


90° 


1 


2.226 


0.278 




18, 19, 


20, 21 












-0.968 


-0.105 


146 + 8 
















-0.709 


-0.507 



















-0.696 


0.439 









SDSS J1004+4112 


0.680 


1.734 


90° 





3.943 


-8.868 




5, 22, 


23, 24 












-8.412 


-0.86 


821.6 + 2.1 
















-7.095 


-4.389 


40.6 + 1.8 
















1.303 


5.327 









WFI J2033-4723 


0.661 


1.660 


0° 





-1.439 


-0.311 




25, 26, 27 












0.756 


0.949 


35.5 + 1.4 
















0.044 


1.068 



















0.674 


-0.5891 


27-l^i 






PKS 1830-211 


0.885 


2.507 







-0.498 


0.456 




28, 29, 


30, 31 












0.151 


-0.268 


26^ 






B1600+434 


0.410 


1.590 


90° 


1 


0.610 


0.814 




5, 32, 


33, 34 












-0.110 


-0.369 


51 + 4 






B0218+357 


0.685 


0.944 







0.250 


-0.119 




35, 36, 


37, 38 












-0.052 


0.007 


lo.ill-j 
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-PG 1115+080 




"SDSSJ1650+4251 






"SBS 1520+530 








SBS 0909+532 





Fig. 1. — Ensemble of 18 average mass maps of the lenses. The larger tick marks in each panel correspond to 
l'/O. Red and cyan dots represent the positions of the images and the sources respectively. The contours are 
in logarithmic steps, with critical density corresponding to the third contour from the outside. The systems 
are grouped according their morphology described in the Section 5. 
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([Oguril 120071 ) and 10 lenses ( Sana et al.l l2QQ6h . 
The Ho distribution of the samples of 10 and 
15 lenses are slightly different than the origin al 
results ofjOguril ([20071 ) and ISaha et all ([20Q6f lEl. 
Oguril (|2007[ ) used 16 lensed quasar systems (40 
image pairs) to constrain the Hubble constant. For 
each image pair, he computed the likelihood as a 
function of the Hubble constant. He then com- 
puted the effective \ 2 by summing up the loga- 
rithm of the likelihoods. The first summation runs 
over lens systems, whereas the second summation 
runs over image pairs for each lens system. He 
derived the best-fit value and its error of Hubble 
constant in the standard way using a goodness-of- 
fit parameter. 



Fig. 2. — Histogram of the ensembles of Ho val- 
ues estimated from 18 lensing systems. We have 
assumed Q\ = 0.7 and Vt rn = 0.3. 



The difference with Oguril (|2007[ ) is mainly due 
to use of different modeling method and a different 
statistical and modeling a p proac h to obtain Ho. 
In the case of ISaha et al.l ([20061 ) the method is 
identical and the difference comes from the use of 
other/newer data, and the use of different rules for 
constraini ng shear, adding s econdary lenses, etc. 
Using the ISaha et all (|2006f ) sample of 10 lenses 
we have o btained Hp = 63^5 km s _1 Mpc _1 and 



using the Oguril ([20071 ) sample of 16 lenses (minus 
B1422+231), we got H = 66^ km s^Mpc" 1 . 



Oguril (120071) ob tained H 



v.,.™. - 70 ± 6 km s _1 Mpc" 

Saha et al.f(l2006h got H = 72± 7 11 km s^Mpc" 1 



Fig. 3. — Comparison of the three histograms of 
the ensembles of Ho values. Histogram plotted 
with the solid line represents our results for 18 
systems; t he dashed line re presents the results of 
10 lenses (I Saha et al.ll2006l) and the dotted line is 
a result of 16 lenses gurill2007l ). All calculations 
are done for a flat Universe with £7 a = 0-7 and 

n m = 0.3. 



5. 'Elliptical' sample 

The strongest degeneracy in lens modeling is 
the so-called mass-sheat degeneracy between time 
delays and the steepness of the mass profile. This 
degeneracy causes that without change in the im- 
ages position, we can change the steepness of the 
mass profile, hence, the resulting Ho. Thus, if 
the steepness and At are known, Ho is well con- 
strained, although the time delays are also influ- 
enced by other more complicated degeneracies in- 
volving details of the shape of the lens but t hese 
effects are secondary ([Saha fc Williams! 120061 ). 

In our sample of 18 systems we have galaxies 
with a variety of steepness. Without detailed ob- 
servations the profile steepness of each lens is not 
known, except possibly for elliptical galaxies. 

Several studies have shown that elliptical galax- 
ies may be consid ered as approximately is other- 
mal p oc r~f ( Koqpmans et al.l 2006 : Oguril 



A p . v 

2007UKoopmans et al]l2009l iGerhard et al.ll200ll ). 
Hence, by selecting from our 18 system only ellip- 
tical galaxies we expect to get a uniform sample 
with known slopes of mass profiles. 

We have selected five systems that are rela- 
tively isolated, ellip tical galaxies: P G 1115+080 



([impev et al.l Il998f ). HE 2149-275 ([Lopez et al 
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Table 2 — Continued 



System 






P.A. 


Point mass a 


x 6 


y b 


Ar c 


Reference 












(") 


H 


(days) 




B1608+656 


0.630 


1.394 


- 


1 


-1.155 


-0.896 




39, 40, 41, 42 












-0.419 


1.059 


31.5 + 1.5 














-1.165 


0.610 

















0.714 


-0.197 


41 ± 1.5 




HE 1104-181 


0.729 


2.319 




2 


-1.936 


-0.832 




43, 44, 45 












0.965 


0.5 


152.21*;° 




SBS 1520+530 


0.761 


1.855 




1 


1.130 


0.387 




46, 47, 48, 49 












-0.296 


-0.265 


130 + 3 




FBQ J0951+263 


0.24 


1.246 







0.750 


0.459 




50, 51 












-0.142 


-0.169 


16 + 2 




SBS 0909+532 


0.830 


1.376 







0.572 


0.494 




5, 52, 53 












-0.415 


-0.004 


4511, 




B1422+231 


0.339 


3.62 


10° 





1.014 


-0.168 




54, 55, 3 












0.291 


0.900 


7.6 + 2.5 














0.680 


0.580 


1.5 + 1.4 














-0.271 


-0.222 








a We have included all the galaxies that might have contribution in the lensing, we add them whenever 
one or more galaxies are visible in the field and when their redshift is similar to the main lens. 

b The positions of the QSO images are calculated with respect to the position of the main galaxy Gl. 
The images are in arrival-time order. 

C A11 measured time delays are listed expect those for which error bars are large and, therefore, the 
detections are marginal. 



References. — (1) IMorgan et all (120081): (2) iBarkanal (Il997h: (3) iTonrv 



(ll9S0|); (5) CASTLES ; 
IOguriet~aTl (12 005) 



113), 



(10) _ 
Vuissoz et all (12 007) 



Falco 



(20) 



Kneibetal 



2003); (24) 



Mor gan et al 



Lidman et al 



(6) iBurud et all (l2002al): (7 ) IWisotzki et al.1 (Il996h : 

Kochanek et al.l (120061) : (11) IMorgan et al.1 (I2005D: 



(1998) 



(4) IWevmann et al 
Paraficz et_al.l (120091): 



(12) 



(14) iBernstein fc Fischer (Il999h : (15) lOscoz et al 



et al.1 (Il997l ): (17) IClaeskens et all I2OO6I): (18) iBumd^al 



(120001): (21) iBade et al.l (Il997l) 



Morgan et ahl 
j200ll)L_a 6 ) 



(Il998l): (lgljHiorth et alJ (l2002l): 
(22) iFohlmeister et al.1 (|2008l): (23 ) llnada et al.1 



Williams fc Sahal (120041): (25) I Vuissoz et al.l (l2008|): ( 26) lEigenbrod et al.1 |2006l): (27) 
(l2QQ4h : (28) lMevlan et al l (120051 ): (29) lLovell et al.l ril998l) ; (30)1chengalur et al.l(| l99 9j); (31 ) 
(Il999l) : (32) IBurud et al.l (l2000l) : (33) Llaunsen fc Hiorthl (Il997l): ( 34) I Jackson et al.1 (Il995l) : 



(35) IWucknitz et all (|2004l): (36) ICohen et al.l (l2000l ): ( 37) iBrowne et al.1 (Il993l) : ICari 
Cohen et al.1 (120031): (39) iK oom nans et al.l (|2003l): (40) l^s snacht et al. (2001 ); (41 ) 



(42) | Fassnacht et al 



li et al.1 dl993^; (38 ) 



Mvers et al.l (Il995l) 



- (Il996l): (4 3) IPoindexter et al l (120071) : (44 ) iLidman et al.l (120001) : ( 45) IWisotzki et al 
( 1993h : (46) lFaure et al.1 (l2002h ^(47) lBurud et al.l (l2002bl )^(48) lAuger et al.1 (l2008h ^49) IChavushvan et al 
( 1997h : (50) [fakobsson et al.1 imm ^l) Iwhite et al.l (I2000l): (52) Ullan et al.l ([20051 ): (53) ILubin et al 
( 2000h : (54) iRavchaudhurv et alT(j2003[ ): (55) IPatnaik fc Narasimhal ([20011 ): 



7 



19981) . SDSS J 1206+4332 (IParaficz et all l2QQ9h. 
HE 0435-1223 (Kocha nek et al.ll2QQ6|) and SDSS 



13 lenses + 5 e 



J1650+4251 (IVuissoz et all l2Q07h Morgan et al. 

Ell). 

Five other lenses were also found to be el- 
liptical but a cluster or group to which they 
belong has strong influence on the lensing sys- 
tem, thus it is difficult to model them: QSO 
0957 +561 (|Oscoz et af]ll997t iBernstein Fischer 
19991 ): RX J 1131-1231 (ClaeskenselaDl2QQ6|); RX 
J0911+055 dKneib et al.l l2QQ0h. SDSS J1004+4112 



I Sharon et al.l2QQ5h and WFI J2033-4723 (jEigenbrod et al.F ig. 4.— Histogram of the distribution of H val- 
20061 ). ues for the selected sample of 5 elliptical galax- 

ies having constrained steepness of mass profiles 
in the range 1.8-2.2 (solid line) and the selected 
elliptical sample combined with the rest of the sys- 
tems (dashed line). All calculations are done for a 
flat Universe with Q\ = 0.7 and ^ m = 0.3. 



Three of the lenses are most probably spi- 
ral gala xies: PKS 1830- 211 is a face-on spiral 
galaxy ([Winn et al.l l2002f ). B1 600+434 is a spi- 
ral g alaxy with a companion ([jaunsen fc Hjorthl 



19971) and B0218+357 is an isolated spiral galaxy 
(|Koopmans fc CLASSll200lh . 



The remaining five systems were not included 
into the 'elliptical' sample due to various other 
issues: B1608+656 has two lensing galaxies in- 
side the Einstein ring; HE 1104+180 has little 
starlight, suggesting: a dark matter dominate d lens 
([Poindexter et al.l 120071 : IVuissoz et al.ll2008l ): SBS 
1520+530 has a steeper than isothermal slope, 
probably due to mergers ([Auger et al.ll2008l ): FBQ 
0951+ 263 is a complicate d system which is hard to 
model (IPeng et alJl2QQ6h: SB S 0909+532 is prob- 



ably early- type (jLehar et al.l l2QQoh but according 



Motta et al.l ((2002) is not very typical due to lots 



of dust. 

Using the sample of five elliptical galaxies and 
constraining the steepness of their mass profiles 



to be dr, 



1.8 and a n 



2.2 we have run 



the PixeLens simultaneous modeling. The sample 
of 5 lensing systems gives us a Hubble constant 
estimation Ho = 79^ km s _1 Mpc _1 at 68% con- 
fidence and H = 79±| km s^Mpc -1 at 90% con- 
fidence. We have also combined the 5 constrained 
systems with the rest of the unconstrained sample 
to perform simultaneous modeling and obtained 
very well determined Hubble constant Ho = 76 + 3 
km s _1 Mpc _1 at 68% confidence and Ho = 76 + 5 
km s _1 Mpc _1 at 90% confidence. The results are 
presented in Figure 4. 



6. Conclusions 

Non-parametric modeling using PixeLens was 
applied to an ensemble of 18 lenses to determine 
a new value for the Hubble constant. We have 
obtained Ho = 66^4 km s _1 Mpc _1 for a flat Uni- 
verse with ft\ = 0.7, Vt rn =0.3. We have also com- 
pared our results with the two p revious attempts 
to estimate Ho from time delays (|Saha et al.ll2006l : 
l Oguril ioog) (Fig. 3). 

Our additional result was based on studies of a 
selected sample of five lensing galaxies that have 
mass profiles close to p oc r -2 . The Hubble con- 
stant recovered using the selected sample of ellipti- 
cal galaxies combined with the rest of the systems 
is H = 76^ km s^Mpc" 1 . 

The gravitational lensing method that con- 
strains Ho has difficulties due to a couple of de- 
generacies between mass and time delay. The 
major degeneracy, the mass-sheet degeneracy, as 
we have shown, can be addressed by a careful 
choice of galaxies and the others partially by a 
combined, pixelated analysis of a large sample of 
lenses. Pixelated lens modeling provides insight 
into the structure of galaxies and the distribu- 
tion of dark matter which together with precise 
measurements of time delays gives a reliable cos- 
mological method. 

Lensing can already determine the Hubble con- 
stant approaching the accuracy level of other lead- 
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ing measurements. Nevertheless, still more obser- 
vation are needed. Future data with precise time 
delay measurements and better lens models will 
give even better constraints on Ho, perhaps turn- 
ing lensing into a very competitive method. 

The Dark Cosmology Centre is funded by the 
DNRF. DP thanks Prasenjit Saha for very exten- 
sive and patient help with PixeLens. 
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